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ABSTRACT 

In this article, a model of random hermitian matrices is considered, in which the mea- 
sure exp(— S) contains a general ?7(iV)-invariant potential and an external source term: 
S = N tr(V(M) + MA). The generalization of known determinant formulae leads to com- 
pact expressions for the correlation functions of the energy levels. These expressions, exact at 
finite N, are potentially useful for asymptotic analysis. 
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1. Introduction. 



As first suggested by Wigner [1] , random matrices can be used to simulate Hamiltonians of 
complex systems. In this approach, one would like to characterize the structure of the energy 
levels, which are represented by the eigenvalues of a large matrix that can be assumed random. 
It is now known that many statistical properties of spectra of true physical systems are indeed 
well described by those of random matrices (cf [9] for a review): it is therefore important 
to understand how much these spectral properties depend on the particular matrix ensemble 
chosen, i.e. determine universality classes of matrix ensembles. 

We shall consider here ensembles of hermitian matrices only, which correspond to sys- 
tems without time-reversal invariance. Let us introduce the probability distributions p n of 
the eigenvalues: If M is a random hermitian N x N matrix, we define p n (Ao, . . . , A n _i) to 
be the density of probability that M has (Ao, • • • , A n _i) among its N eigenvalues, with the 
normalization convention that: J Yi^o d\ Pn(Ao, • • • , A n _i) = 1. Following [2], we also define 
R n = (J^J n \\ Pn, a different normalization which makes the connection with the correlation 



since these quantities only differ by 5 functions for coinciding eigenvalues. 

In the case of a U (iV)-invariant measure, consisting of a simple potential term, which we 
briefly review in section 2, the theory of orthogonal polynomials allows to write down exact 
expressions at finite iV for these functions, in terms of a single kernel K(X,/i) [2]: 



The study of the distribution of eigenvalues then boils down to the analysis of this kernel; 
in particular the asymptotics of K as N — > oo allow to compute asymptotics of correlation 
functions and find the different "universal" properties that arise in this limit. Therefore it 
seems quite interesting to find similar formulae for more general measures. In fact, in recent 
papers [12], Brezin and Hikami have shown that formula (1.2) can be generalized to the gaussian 



functions 




for distinct Aj 




-Rn(Ao, Ai, . . . , A n _i) = det(K(Xi, Xj))ij=o... n -i- 



(1.2) 
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ensemble with an external field. More precisely, for the (unnormalized) measure 

N ■ ,1 AT*- 1 / a\ jN 2 



exp ^ — — trM + NtiMAj d M, (1.3) 
they introduced the kernel 

K(X u) ~ — I ' —i—T\ ( lt / N ~ ai \ I ^-f^N-itX+Nvv (1 4) 

1 ' W N J 2rr J 2m f-} Q \ v - a t Jit/N-v 1 J 

(the ~ is here to distinguish this kernel from a slightly different one that will be introduced 
later) where the ai are the eigenvalues of the hermitian matrix A and the contour integral 
encircles these eigenvalues. (1.2) then holds with K instead of K. 

The aim of this paper is to define a kernel K such that (1.2) still holds in the more general 
case of an arbitrary potential with an external source term, which is the subject of section 
3. As an example we consider in section 4 the gaussian ensemble with an external field and 
reproduce the kernel K obtained in [12]. 

2. The U (N) -invariant case. 

Let us consider an ensemble of random hermitian N x N matrices with the measure 

Z" 1 exp(-iVtr\/(M))rf iv2 M (2.1) 

where V is a polynomial and Z the partition function. A classical result [2] expresses the 
distribution law p n of n eigenvalues (1 < n < N) of M in terms of the kernel 

N-l 

K(X,fi)= ^F k (\)F k (iM). (2.2) 

Here Fj is the orthonormal function associated to the usual orthogonal polynomial -Pi (A) = 
A* + ---: 



F(X) = h7 1 / 2 P t (X)e W (-^V(X)^ 
J dAexp(-W(A))Pi(A)Pj(A) = h^. 



(2.3) 



(see [11] for a review of orthogonal polynomials in matrix models). Let us briefly rederive 
this result in a manner that naturally generalizes. As the measure (2.1) only depends on the 
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eigenvalues of M, the integration over the angular variables is trivial and one finds: 

PN (\ , Ai, . . . , Aiv-i) = Z- 1 A\X i ) exp (-N ^ V(Xi) j . (2.4) 

The Van der Monde determinant A(A^) = det(Aj : ')j i j = o...jv-i can be rewritten in terms of the 
orthogonal polynomials: 

/ N-l 

Piv(Ao, Ai, . . . , Ajv-i) = Z^ 1 det(P A; (Aj)) i;A;= o...jv-i det(P&(Aj)) J>fc=0 ...jv-i exp I -V V(Xi 

V i=0 

(2.5) 

One can now easily compute Z = N\ n^o 1 ^* by integrating over all Aj. Combining the two 
determinants, we finally obtain: 

Piv(Ao, Ai, . . . , Ajv-i) = det(K(Aj, \j))ij=o...N-i- (2.6) 

The kernel if has the following properties: 

K(\,ri = K(»,\) 



[K*K\(\,p) = J dpK(X,p)K(p,p) = K(X,p) 



(2.7) 



i.e. it is the orthogonal projector on the subspace spanned by the F^, < k < N — 1. Using 
the property K -k K = K and noting that 

Pn(Ao, Ai, . . . , A n _i) = J d\ n p n+ i(\ , Ai, . . . , A n ), (2.8) 

one can then show inductively that 

(N — n)\ 

p n (A , Ai, . . . , A n _i) = — — — det(K(\i, Xj))ij=o... n -i- (2.9) 
for any n < N . This is equivalent to formula (1.2). 



3. Generalization to the case of an external field. 

We shall now see that in the case of a general measure with an external field, (1.2) still 
holds; a simple expression for a kernel K will be derived. Let us indeed consider the measure: 

Z- 1 exp(-NtrV(M) + NtrMA)d N2 M (3.1) 

where V is an arbitrary polynomial, and A = diag(ao, . . . , a/v-i) can be assumed diagonal. 

Particular matrix models of this type and their large N study appear in many papers [5,7]. 
Here we shall go beyond the 1/iV-expansion and write exact expressions at finite N. 

One diagonalizes M: if M — QAQ} where A = diag(Ao, . . . , Ajv-i), the integral over Q is 
the usual Itzykson-Zuber integral [4] on the unitary group and we find: 

p N (Xo, A b . . . , A*_0 = Z-^AM **^™^ exp (-N g V()^J . (3.2) 

We replace as usual powers of A in the Van der Monde with the orthogonal polynomials P k (X) 
of the measure exp — NV(X)dX. Z can now be computed: 

N-l N-l 

Z = m AT~) / 11 d\ i det(P k (\ i ))expNj2(-V(^) + ^) 

[a ° *=0 i=0 (3.3) 

= det ( [ d\P k (X) exp N(-V{\) + ai X)) 

A ( a U \J Jk,l=0...N-l 

Inserting (3.3) into (3.2) yields 

p N (\ , A b . . . , Aiv-i) = —r TTTTTTFrn ^ T77TVT TvT exp _iV 2^ y ( Ai ) 

JV!det(J dXP k (X)expN(-V(X) + a,A)) fc/=0 jV _ 1 ^ ^ 

(3.4) 

This formula has a remarkably simple structure, which can be made more explicit by intro- 
ducing F fc (A) = h k 1/2 P k (X) exp-f F(A) as before, and Gj(A) = exp [iVa/A - f V(X)]: 

/ \ \ \ \ 1 det(F fc (Aj)) i;fc=0 ...iv-i det(G{(Aj)j/ = o...iv-i / Q _v 

' A) """ i «- ,I = S det(/dAF t (A)G,(A)W..*-i ' (3 ' 5) 

The matrix (J dXGi(X)F k (X))i jk= o ...jv-i possesses an inverse, which we denote by a^; putting 



together the three determinants we finally obtain: 



PivOo, Ai, . . . , Ajv-i) = det(K(\i, Xj))ij= ...N-i (3.6) 

where 

N-l 

K(X,fi) = F k(^hklGi(fi). (3.7) 

k,l=0 

The kernel X satisfies the property: 



N-l r „ 

[K*K\{\,p)= a kiF k {\) dfiG^F, 

h 1.1 1 li— n \-J 



k,k',l,l'=0 
= K(\,p). 



®k'i>Gi'(p) 



(3.8) 



Thus, one can follow the same line of reasoning as in the ?7(./V)-invariant case to obtain the 
determinant formulae 



(TV — n)\ 

Pn(Ao, Ai, . . . , A n _i) = — — — det (K(Xi, Xj))i,j=o... n -i 



(3.9) 



for any n < N . 

The kernel K is of the form 



N-l 

K(X,p) = J2 F k( X ) P M (3.10) 

fc=0 

with -Ffe(/i) = J2i a klGi(fj,); but F k ^ F^ and therefore K is not symmetric. Further analysis 
of K resides in the observation that as a — > 0, 

/]\rk , . 

dAP fc (A)expiV(-\/(A) + aA) = /i fc — a fc + O (a fc+1 J (3.11) 

due to the orthogonality of the P^. The first consequence is that it is only in the limit A — > 
that the G/ become linear combinations of the F^, < k < N — 1, ensuring that Ffc — > as 
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expected. The second consequence is that the quantity 



N-i r 

a~ N Y[ (a ~ ai) / dXP k (X) exp N(-V(X) + aX) (3.12) 
1=0 J 



is well-defined for any k > N; we can define the inverse Laplace transform 



N 



$ fc = exp -V(X) d- N H (d - a,) [P fc (A) exp -NV(X)] , (3.13) 



(where d = \/N d/dX) which is a regular function with fast decrease at infinity. One can see 
that we now have a complete set of eigenvectors of K: 

I K(\,n)F k (n)dn = F k (X) Vk<N 

' (3.14) 
K(X,fi)$ k (fi)dfi = \/k>N. 



I 



Again, it is only when A = that & k = F k . 

We conclude that K is a non-orthogonal projector on the space spanned by the F k) < 
k < N- 1. 



4. Example: the gaussian ensemble. 

In the case where V(M) = \M 2 , the kernel K defined in eq. (3.7) should be related to the 
kernel K of [12]. This is indeed what we find. 

The orthogonal polynomials Pfe(A) are Hermite polynomials, their normalization is h k = 
a/27T / 'Nk\/N k , and one can compute explicitly the integral 

j dXP k {X)e N{ -^ +a ^ = y^jNa k e^ a2 (4.1) 

in accordance with eq. (3.11). One can also directly calculate Z by changing variables from 
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M to M — A in the gaussian matrix integral: this immediately gives 



Z~ ew -J2al (4.2) 
1=0 

up to A-independent factors. We deduce the explicit form of the kernel K: 

rw N - 1 

K{\ ^ = J^H P fc (A)^e^^-K A2 +^)-^?). (4.3) 

71 k,l=0 

where /3m = a^i \f2ix /N exp Y a f is the inverse of the Van der Monde matrix (a/ fc ). By definition 
one has: 

to v# ai - av 

which allows to express fikl m terms of symmetric functions of the (a/')zyi- We now use an 
inverse formula of eq. (4.1): 

+ioo 

P k (X)= [ -^aWW-^+W (4.5) 
J i\/2nN 



where the integration is along the imaginary axis. Combining eq. (4.3), (4.4), and (4.5) we 
obtain: 

N-i , 

K(\,v) = J2 / ^\{^^e N ^ a2 - aX+ ^ 2 --y-^ +a ^. (4.6) 
'=0 -ioo 2m 'Vi ai ai ' 

This can in turn be represented by a contour integral in the complex plane which picks up 
poles at v — af. 

+ico N—l 

- [ £ f£ n fe) -L-^^k-w-^). (4.7, 



2m J 2m - LJ - \v — ai I a — v 

—o v 7 

The redefinition 

K(X, fj) = e-^ x2 K(X, /i)e + ^ 2 , (4.8) 
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gives 

+ioo N—l 
-ioo l ~° 

which coincides with eq. (1.4) if one sets a = it/N. Notice that transformation (4.8) does not 
affect the value of determinants of the type (1.2). 



5. Conclusion. 

We have derived determinant formulae for the correlation functions of eigenvalues in terms 
of a kernel K which has very simple algebraic properties: it is a projector on a space of 
dimension N; its eigenvectors are known, even though they are not as simple as in the U (-AO- 
invariant case, since K is no longer symmetric. 

For the simple gaussian ensemble, the asymptotic behavior of the kernel K{\, ji) when 
A — pi is of order l/N, N — > oo, has been known for a long time [2,3]: 

K(\,(j,)~—, s~N(\-n). (5.1) 
s 

It is remarkable that in the general £7(iV)-invariant case [8] as well as in the gaussian case 
with an external source [12], this behavior remains identical; combined with the determinant 
formulae, it implies universal properties for the correlation functions, as well as for the level 
spacing distribution. The latter quantity, first introduced by Wigner [1], is empirically known 
to be universal for a wide range of models. Here, it can be expressed in terms of K using 
the determinant formulae, and the relation takes a particularly simple form in the domain 
of universality, that is when one considers intervals of order l/N: if p(s) is the level spacing 
distribution, the spacing variable 9 being defined by s = N9, one has 

p(s) = -^det(l-k) (5.2) 

where K is the asymptotic form of K in the interval [—9/2,9/2] as N — > oo, 9 — > 0, s fixed 
(cf eq. (5.1)), and det is the usual Fredholm determinant. Thus, the universality of K implies 
the short distance universality of the level spacing. As opposed to the universal behavior of 
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the correlation function p2(X,p) when A — /i is of order 1, which can be obtained by various 
standard large N techniques [6,8,10], the short distance universal behavior is more difficult to 
find, and its generalization to an arbitrary potential with an external source (using the newly 
found kernel K) is currently under study. 
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